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Fast Finite Field Hartley Transforms Based on Hadamard Decomposition 
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Abstract 

A new transform over finite fields, the finite field Hartley transform (FFHT), was recently introduced and a number of promising 
applications on the design of efficient multiple access systems and multilevel spread spectrum sequences were proposed. The FFHT 
exhibits interesting symmetries, which are exploited to derive tailored fast transform algorithms. The proposed fast algorithms are 
based on successive decompositions of the FFHT by means of Hadamard-Walsh transforms (HWT). The introduced decompositions 
meet the lower bound on the multiplicative complexity for all the cases investigated. The complexity of the new algorithms is 
compared with that of traditional algorithms. 
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1 Introduction 

Discrete transforms defined over finite fields, such as the finite 
field Fourier transform (FFFT), pivotal tools in coding theory |2| 
and signal processing m. Another interesting example is the 
finite field Hartley transform (FFHT), a self-inverse transform 
(involution operator) introduced in |[3l [4l |5). Recent promising 
applications of discrete transforms concern the use of the FFHT 
to design digital multiplex systems, efficient multiple access sys¬ 
tems {6} and multilevel spread spectrum sequences (7). A deci¬ 
sive factor for applications of discrete transforms has been the 
existence of the so-called fast transforms (FT) for computing it. 
Since the FFHT is a more symmetrical version of discrete trans¬ 
form, in this paper this symmetry is exploited so as to derive 
new FTs that require less operations. These FTs, derived for 
short blocklengths (N < 24), are based on successive decompo¬ 
sitions in a similar way as the multilayer Hadamard decompo¬ 
sition employed to compute the discrete Hartley Transform 
(DHT) (9). This new approach, which is based on decompo¬ 
sition of the FFHT by means of Hadamard-Walsh transforms 
(HWT), meets the lower bound on the multiplicative complexity 
of a discrete Fourier transform (DFT) IflOl . Each HWT imple¬ 
ments pre-additions and post-additions. These schemes are easy 
to implement using digital signal processors (DSP) or low-cost 


high-speed dedicated hardware. The complexity of these new 
FTs is compared with that of traditional methods, such as the 
Cooley-Tukey radix-2, split radix, Winograd, and Rader-Brenner 
algorithms, which were adapted to compute the FFHT ifTH . 

2 The Finite Field Hartley Transform 

Finite field Hartley transforms are based on a trigonometry over 
Galois Fields GF(q), q = p r , p = 3 (mod 4), so that ( p — l) 1 / 2 ^ 
GF(q). The set G(q) of Gaussian integers over GF(q) plays an 
important role in this analysis. This set defines a structure GI(q), 
which is isomorphic to GF(q 2 ) |f3l . 

Definition 1 Let C be an element of Gliq) with multiplicative 
order N, where q = p r . The trigonometric functions sine, co¬ 
sine, and cas ( cosine-and-sine or Flartley kernel) are defined by, 
respectively: 

n - r-i 

sin («) = ^^-, 

2 J 

n + C~‘ 

cos(Z) = 

cas(i) = sin(i) +cos(i), 
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for i = 0,1,... ,N — 1. 

Definition 2 Let v = {vo, vi, ..., vjv-i } be a vector of GF(q)- 
valued components, q = p r . The finite field Hartley transform 
(FFHT) is the vector V = {VojVi ,• • ■,V)v-i}> with components 
14 € GI(q m ) given by 14 = v < - cas (ik), where Q is an ele¬ 
ment of multiplicative order N over GI(q' n ). 

The inverse FFHT is given by the following theorem. 

Theorem 1 The vector v = {vo,vi,..., vjv_ 1 } can be derived 
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Table 1: Minimal Multiplicative Complexity achievable for the 
/V-point DFT 


N 4 8 12 16 

H(DFT(N)) 0 2 4 10 


It follows that: 


rvoi 


rO 1 0 11 


Ao(lf 

Vi 


6 0 10 


Sl(l) 

v 2 


0 6 0 1 


s 2 (i) 

LvsJ 


L l o l oJ 


s 3 ( i) 


The number of additions is reduced from 12 to 8 (4 pre-additions 
and 4 post-additions). 


from its FFHT according to: 


N 


1 

(mod p) 


N -1 

£ Vfc-cas (ik), 
k =0 


3.2 Computing the 6-point FFHT 

Let v <—> V be an FFHT transform pair over G/(7). Considering 
£ = 3, the FFHT can be computed by 


for i = 0,1,... ,N — 1. 

3 Hadamard Decomposition of the FFHT 

The Hadamard decomposition was employed in (8) as a tool to 
compute the discrete Hartley transform. This approach allows 
the minimization of the multiplicative complexity of the DHT for 
some blocklengths. Since all the properties and symmetries of 
the DHT are also observed for the FFHT, the application of this 
algorithm to finite fields should be expected. The minimal multi¬ 
plicative complexity of a DFT with blocklength N —denoted by 
p(DFT(N ))—can be calculated by converting the DFT in a set 
of cyclic convolutions. A lower bound on p(DFT{N)) is pre¬ 
sented in ifTOl . Table Q] shows a few values of n(DFT(N)) for 
short blocklengths. 

Considering with finite held transforms, the following com¬ 
ments are worthwhile: 


fVoi 


Vi 


v 2 


v 3 


v 4 


Lv 5 J 



it ill l 

1 4 +j 3+j 6 3+6 ;'4+6 ;' 
1 3+j 3+6 j 1 3+j 3+6 j 
16 16 1 6 

I 3+6 j 3+j 1 3+6; 3+j 
1 4+6; 3+6; 6 3+;' 4+;' 



rvoi 


Vl 


V2 


v 3 


v 4 


i v 5J 


Observing the symmetries, a first column combination can be 
made. Let 


So(l)=(v4-vi),Si(l)=(v4+v 1 ),S 2 (l)=(v 5 -V2), 

S3(1)=(v'5+''2)A 4 (1)=(vo-V3),S 5 (1)=(vo+V3). 


Therefore, 


fVoi 


r o 

1 

0 

1 0 1-1 


rSo(i)i 

Vl 


3+6; 

0 

4+6; 

0 1 0 


Sl(l) 

v 2 


0 

3+j 

0 

3+6; 0 1 


*(1) 

v 3 


1 

0 

6 

0 1 0 


S 3 (l) 

v 4 


0 

3+6 j 

0 

3+; 0 1 


S 4 (t) 

Lv 5 J 


L 3+; 

0 

4 +; 

0 1 0. 


L5 5 (1)J 


Going on with this procedure, a second pre-addition layer is de¬ 
rived: 


(i) The minimal multiplicative complexity, ji(FFFT(N)), 
for a FT over the finite held GI(p r ), is the same as 
H(DFT(N)), evaluated over the real held. 

(ii) The relationship between the multiplicative and additive 
complexity over a hnite held strong depends on implemen¬ 
tation. For small p. the total complexity (additive plus mul¬ 
tiplicative) must be taken into account since their difference 
is small. 

New algorithms for computing the FFHT are introduced in the 
next section. 

3 . 1 Computing the 4-point FFHT 

Let v 4—* V be a FFHT transform pair over GI( 7). The FFHT, 
assuming a cas(-) kernel with £ = j, is computed by: 


rvoi 


ri i i n 


pol 

Vl 


116 6 


Vl 

V 2 


16 16 


v 2 

Lv 3 J 


L 1 6 6 1 J 


v 3 


So(2) = S 2 (l) — So(l),Si(2) = S 2 (l) +S 0 (1), 
S 2 (2)=S 3 (1)-S 1 (1),S3(2)=S 3 (1) + S 1 (1), 
S 4 (2) = S 4 (1),S 5 (2) = S 5 (l). 

Finally, 


r^oi 


rO 0 0 10 1-1 


rSo(2)-i 

Vi 


4 6; 0 0 1 0 


5l(2) 

v 2 


0 0 6; 3 0 1 


52 (2) 

v 3 


6 0 0 0 1 0 


5 3 (2) 

V 4 


0 0 ; 3 0 1 


S 4 (2) 

Lv 5 J 


.4 1; 0 0 10. 


L^ s ( 2 )J 


Since there is only one multiplication (by the same factor) in 
columns 1 and 4, there will be two multiplications. The total 
number of additions required to compute a 6-blocklength FFHT 
is 16 (10 pre-additions and 6 post-additions). 

3.3 Computing the 8-point FFHT 

Let v i—> V be an FFHT transform pair over G/(7). Let C = 
2+2j, so the corresponding matrix formulation is, 


Indeed, no multiplication is needed. Observing further symme¬ 
tries, columns can be combined through Hadamard blocks in or¬ 
der to reduce the number of additions. Let 
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lv 7 j 
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11661166 
10646013 
16 16 16 16 
13106460 
1 6 6 1 16 6 1 
1 0 6 3 6 0 1 4- 
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v 4 
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Lv 7 J 


So(l) = (v 3 -vi),Si(l) = (V3 + V1), 
£>(1) = (vo — v 2 ), 5 3 (1) = (vo + v 2 ). 
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Defining a 1st order pre-addition layer: 

So(l) = (v5-Vl),Si(l) = (V5 + V1), 
S 2 (l) = (l ; 6-V2),S 3 (l) = (V6 + V2), 
54(1) = (v 7 — V3), 55(1) = (V7 + V3), 
^6(1) = (v 0 — V 4 ), 57(1) = (V 0 + V 4 ). 

Therefore, 


rv 0 i 

Vi 

v 2 



rSo(i)-| 
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Sl(l) 
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01060601 


S 2 ( 1) 
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s 4 (l) 
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5 5 (1) 

*56(1) 



L<5 7 (1) J 


Defining a 2nd pre-addition layer, 

So( 2 ) = S 0 ( 1 ),S 1 ( 2 )=S 4 ( 1 ), 

S 2 (2) = S5(1)-S 1 (1),S3(2)=S 5 (1) + S 1 (1), 
5 4 (2) = S 6 (l) — S 2 (l), S s (2) = S 6 (l) + S 2 (1), 
S 6 { 2) = S 7 (l)-S 3 (l), S 7 (2) = 5 7 (1) +5 3 (1). 

Consequently, we obtain: 
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<S 6 ( 2 ) 
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Since there is only one multiplication (by the same factor) in 
columns 1 and 2, there are two multiplications. The number 
of additions is 22 (14 pre-additions and 8 post-additions). It is 
worthwhile to remark that the additive complexity is less than 
the one for an 8-DFT calculation by the Winograd algorithm. 

3.4 Computing the 12-point FFHT 

Let v i —> V be an FFHT transform pair over GI(p). As an ex¬ 
ample, let p = 7 and £ = 3 j. Then V = Tv where, 

'l 1 111 111 111 1 ■ 

1 4+63 4 + 6 ; 1 3 + 6 3 4+3 6 3+3 3+3 6 4+3 3 + 6 3 

I 4+63 3+6 3 6 3+3 4+3 1 4+63 3+6 3 6 3 +3 4 +3 

II 661 166 116 6 

1 3+6 3 3 +y 1 3+6 3 3+3 1 3 + 6 ; 3+3 1 3+6 3 3+3 

T _ 1 4+3 4+3 1 3+3 4+63 6 3+63 3+63 6 4+63 3+3 

J ~ 16 161 616 161 6 
1 3+3 4+63 6 3+63 3+63 6 4+63 3+3 1 4+3 4+3 

1 3+3 3+63 1 3+3 3 + 6 3 1 3+3 3+63 1 3+3 3+63 

16 611 661 166 1 
1 4+3 3+3 6 3+63 4+63 1 4+3 3+3 6 3+63 4+63 
.1 3+63 4+3 6 3+3 3+3 6 4+3 3+63 1 4+63 4+67. 

Defining a 1st order pre-addition layer, we obtain: 

So(l) = (v 7 —vi),Si(l) = (V7 + V1), 

S 2 (l) = (V 8 - V 2 ), 53 ( 1 ) = (V 8 +V 2 ), 

S 4 (l) = (v 9 -v 3 ),S 5 (l) = (v 9 + v 3 ), 

^ 6 ( 1 ) = (' ; 10 — l ; 4 )i S 7 (l) = (V10 + V4), 

S 8 (l) = (vn — V5), 5 9 ( 1 ) = (vn +V5), 

■Sio(l) = (i ; 0 — V6), 5 il(l) = (vo + V6). 


Therefore, V = where 
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A second order pre-addition layer can be defined according 

S 0 (2) = S 6 (l) -S 0 (l), S 6 (2) = S 9 (l) —S 3 (l), 

Si(2) = S 6 (l) + So(l),S 7 ( 2 )= 59 (l)+ 5 3 (l), 

S 2 (2) = S 7 (l) -Si(l), S 8 (2) = Sio(l) - S 4 (l), 
S 3 (2)=S7(1) + Si(1),S 9 (2)=Sio(1) + S 4 (1), 
S 4 (2)=S 8 (1)-S 2 (1),S 1 o(2)=S 11 (1)-S5(1) j 
S 5 (2)=S 8 (1)+S 2 (1),S 11 (2)=S 11 (1)+S 5 (1). 

Therefore, V = T^S(2), where 
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Going further, a 3rd order pre-addition layer is defined: 

So(3) =S 4 (2) -Si ( 2 ), S 6 (3) = S 7 (2) — S 3 ( 2 ), 

Si(3) =S 4 (2) +5i(2), 5 7 (3) = S 7 (2) +S 3 (2), 

S 2 (3) =S 5 (2) —5 0 (2), 58(3) = S 8 (2), 

S 3 (3) =^ 5 ( 2 ) +S 0 (2), 5 9 (3) = S 9 (2), 

S 4 (3) =S 6 (2) -S 2 (2), Si 0 (3) = S 10 (2), 
S 5 (3)=S6(2)+S 2 (2),Sh(3)=S 1 i(2). 

Thus V = T^S( 3), where 

"000000 0 10001 " 

63 03000 001000 
00004 3 000010 
060000 000100 
000000 330001 
t(3) _ J 03000 001000 
1 — 000060000010 • 

040 3 00 000100 
000000 67 30001 
00 1 000 001000 
00004 63 000010 
. 0 4 6 / 0 0 0 0 0 0 1 0 0. 

Since there is just one multiplication by the same factor in 
columns 2, 3, 5 and 8, the total number of multiplications is 4. 
The number of additions required to compute the FFHT is 44 (32 
pre-additions and 12 post-additions). The multiplicative com¬ 
plexity reaches the minimum theoretical complexity and again 
the additive complexity is the same as the one obtained for the 
DHT |D. 
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3.5 Computing the 16-point FFHT 

Let v' 4 —> V a FFHT pair over GI(p). Assuming p = 7 and 
£ = 2 + 4/, the corresponding transform is V - T v, where 


liiiiiiiiiiitiii- 

1 2 j 4 2 j 1 6 j 0 j 6 3/ 3 2 j 6 j 0 67 
1410636014106360 
1 2/ 0 5; 6 67 4 6 j 6 5; 0 2; 1 7 3 j 
1166116611661166 
1 6; 3 6; 1 5 j 0 2 j 6 ; 4 j 6 2 j 0 5; 
1064601310646013 
1 6; 0 6; 6 2j 3 2; 6 67 0 j 1 5; 4 5; 
1616161616161616 
1 5/ 4 5; 1 j 0 6 j 6 2 j 3 2j 6 6 j 0 j 
1310646013106460 
1 5 j 0 2; 6 / 4 j 6 2; 0 5; 1 67 3 6; 
1661166116611661 
1 j 3 j 1 2j 0 5j 6 6; 4 67 6 5j 0 2/ 
1063601410636 0 14 
1 6; 0 7 6 5j 3 5./ 6 j 0 67 1 2./ 4 2/. 


Then, V = T^S(2), where 


T = 


Now we consider the first order pre-addition layer according to: 

So(l) = (v 9 — vi), S'i(l) = (v 9 + vi), 

^2(1) = (vio + v 2 ), ^3(1) = (vioT V2); 

5 4 (1) = (on - v 3 ), 5 5 (1) = (on +V3), 

^6(1) = (l'12 - V4), £7(1) = (l'12 + V4), 

^b(I) = ( 1 '13 — V5), s 9 (l) = (l'13 TPs)> 

Sio(l) = ( 1 '14 — 1'6), 5n(l) = (V14 + V6), 

^12(1) = (l'l5 ~ V7), 5"i3 (1) = (V15 +V7), 

^14(1) = (l'0 — P8), 5 i 5 ( 1) = (l'o + Vs). 

Therefore, V = T [ 1 - 1 .S’ (1), where 


7T 1 ) = 


010101010101010 1 - 
5 j 0 3 0 5 j 0 6 0 j 0 0 0 67 0 1 0 
0401 000603060001 
57 0 0 0 27 0 1 0 7 0 3 0 j 010 
0106060101060601 
j 04 0 7 0 6 0 27 0 0 0 57 0 1 0 
000604060001 0301 
67 0 0 0 7 0 1 0 57 0 4 0 57 0 1 0 
0601060106010601 
27 0 3 0 27 0 6 0 67 0 0 0 j 010 
0301 000604060001 
2 j 0 0 0 57 0 1 0 67 0 3 0 67 0 1 0 
0606010106060101 
67 0 4 0 67 0 6 0 57 0 0 0 2 j 0 1 0 
000603060001 0401 
.7 0 0 0 67 0 1 0 2 j 0 4 0 27 0 1 0. 


By defining a 2nd order pre-addition layer, we have: 


So(2)=S 4 (1)-S 0 (1),Si( 2) = 
S 2 (2)=S 9 (1)-S 1 (1),S 3 (2) = 
S 4 (2)=S 2 (1) } S5(2)=S 10 (1), 
S 6 (2)=5 11 (l)-5 3 (1),5 7 (2) = 
S 8 (2)=Si2(l)-S 8 (l),S 9 (2) = 
S 10 (2)=S 13 (1)-S 5 (1 ),Si 1 (2) 
S 1 2(2)=Si 4 (1)-S 6 (1),Si 3 (2) 
S 14 (2)=Si 5 (l)-S 7 (l),Si5(2) 


S 4 (l) + S 0 (l), 
S 9 {l) + S!(l), 

= S n (l)+5 3 (l), 
= Si 2 (1)+S 8 (1), 
= S 13 (1)+S 5 (1), 
= 5i 4 (1)+5 6 (1), 

= Si 5 (1)+S 7 (1). 


p( 2 ) = 


■ 0 0010001 0 
0 57 0 0 3 0 0 0 67 
0 03000600 
27 00003000 
0 00100060 
0 j 0 0 4 0005; 
0 00000100 
7 0000400 
0 0060001 
0 27 0 0 3 0 0 0 
0 0400060 
5 j 0 000300 
0 00600060 
0 67 0 0 4 0 0 0 27 
0 00000100 
.6/ 0 00040 0 0 


ooioooi- 
0001000 
0000010 
j 0 0 0 1 0 0 
0060001 
0001000 
0 3 0 0 0 1 0 
57 0 0 0 1 0 0 
0060001 
0001000 
0000010 
67 0 0 0 1 0 0 
0010001 
0001000 
0400010 
2 j 0 0 0 1 0 0. 


(1) 


The columns do not cope. However, multiplying both the 5 and 
6 columns by 2 G GF(1), they can be combined with columns 
13 and 14, respectively. Defining then a 2nd order pre-addition 
layer (with two multiplications in columns 10 and 11), we have: 

So(3)=S 0 (2),Si(3)=S 1 (2), 

S 2 (3)=S 8 (2),S 3 (3)=S 9 (2), 

S 4 (3)=S 2 (2),S5(3)=S 10 (2), 

S6(3)=Si 1 (2)-S 3 (2),S7(3)=S 11 (2)+S 3 (2), 

Ss(3) = Si 2 (2) —4S 4 (2), S 9 (3) = S 12 (2) +25 4 (2), 

Sio(3) = Si 3 (2) -4S 5 (2),Sh(3) = Si 3 (2) +25 5 (2), 

S 12 (3) = S 14 (2) -S 6 (2), S 13 (3) = S 14 (2) +S 6 (2), 

Su( 3) = Sis(2) -s 7 (2), Sis(3) = S l5 (2)+S 7 (2). 

Finally, V = T^ 3 ^5'(3), where 


T ( 3 ) = 


000 
0 57 67 
000 

2 j 0 0 
0 0 0 
0 7 57 
0 0 0 
7 0 0 
0 0 0 
0 27 j 
0 0 0 
5 j 0 0 
0 0 0 
0 67 2 j 
0 0 0 
67 0 0 


0 0 0 0 1 0 0 0 0 0 0 0 1 - 
0 000010000000 
0 300000001000 
7 000000100000 
0 006000000010 
0 000001000000 
0 030000000100 
57 000000010000 
0 000600000001 
0 000010000000 
0 400000001000 
67 000000100000 
0 001000000010 
0 000001000000 
0 040000000100 
27 00000001000 0. 


( 2 ) 


There is only one multiplication in columns 1, 2, 3, 4, 5 and 
6, besides two pre-multiplications, so the total complexity is 8. 
The number of additions is 50 (40 pre-additions and 16 post¬ 
additions). In this case, the number of multiplications is less than 
10, the minimum expected multiplication complexity (10]. It can 
be concluded that there are two trivial multiplications. It is not 
simple to identify which are the trivial multiplications from the 
observation of matrices over GI(1). Carrying on the same anal¬ 
ysis over another finite field, the same combination of columns 
was observed, i.e. the approach does not depend on the finite 
field but on the length. Let v 4 —> V be an FFHT pair over GI(p). 
Considering now p = 31, £ = 7+13 j, the transform matrix T 
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will be. 


rl 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -\ 

1 20 0 11 30 6 23 6 30 11 0 20 1 25 8 25 

1 0 30 23 30 0 1 8 I 0 30 23 30 0 1 8 

1 11 23 11 1 25 0 6 30 20 8 20 30 6 0 25 

1 30 30 1 1 30 30 1 I 30 30 1 1 30 30 I 

1 6 0 25 30 11 8 11 30 25 0 6 1 20 23 20 

1 23 1 0 30 8 30 0 1 23 1 0 30 8 30 0 

T _ 1 6 8 6 1 11 0 20 30 25 23 25 30 20 0 11 

1 ~ 1 30 1 30 1 30 1 30 I 30 1 30 1 30 1 30 ' 

1 11 0 20 30 25 23 25 30 20 0 11 1 6 8 6 

1 0 20 8 30 0 1 23 I 0 30 8 30 0 1 23 

1 20 23 20 1 6 0 25 30 11 8 11 30 25 0 6 

1 1 30 30 1 1 30 30 I 1 30 30 1 1 30 30 

1 25 0 6 30 20 8 20 30 6 0 25 1 11 23 11 

1 8 1 0 30 23 30 0 1 8 1 0 30 23 30 0 

L1 25 8 25 1 20 0 11 30 6 23 6 30 11 0 20-1 

Considering the first order pre-addition layer (Equation Q} and 
the second order pre-addition layer (Equation^, V = r f' 2, S(2), 
where 

r 0 0 0 1 000 1 000 1 0 0 0 In 
20 00008000 25 000100 
000000 1000800010 
0 20 0 0 8 0 0 0 25 0 0 0 1 0 00 

0 0 0 30 0 0 0 30 0 0 0 10001 

6 0 0 0 0 23 0 0 0 20 0 00100 

008000 30 000000010 
T ( 2 ) _ 0 25 0 0 23 0 0 0110 0 01 0 00 

1 ~ 000 30 0001000 30 0001 ' 

11 000080006000100 
000000 1 000 23 00010 
01100800060001000 
000 1 000 30 000 30 0001 
25 0 0 0 0 23 0 0 0110 00100 
00 23 000 30 000000 0 10 
L 0 6 0 0 23 0 0 0 20 0 0 0 1 0 0 0-1 

Again, some columns do not cope. But, multiplying both 
columns 5 and 6 by 4 € GF (31), they can be combined with 
columns 14 and 13, respectively. The same occurs with columns 
9 and 10, which combine with columns 2 and 1, respectively. A 
third layer of pre-additions, including four pre-multiplications in 
columns 1, 2, 5, and 6, is given by: 


So(3) = S 7 (2) - S 3 (2), Si (3) = S 7 (2) + S 7 (2), 

S 2 (3) = S 8 (2) -7Si(2), S 3 (3) = S 7 (2) + 8Si(2), 

S 4 (3) = Sg(2)-7S 0 (2), S s (3) = Sg(2) + 7S 0 (2), 

S 6 (3) = S 2 (2),S 7 (3)=Si 0 (2), 

S 8 (3) = S 12 (2)-4S 4 (2),S 9 (3) = S n (2) + 4S 4 (2), 

Sio( 3) = S 13 (2) -4S 5 (2), Sh(3) = S 13 (2) +4S 5 (2), 
S 12 (3) = S 14 (2) - S 6 (2), S 13 (3) = S 14 (2) + S 6 (2), 

S 14 (3) = Si 5 (2) - S n (2), S 15 (3) = S 15 (2) +S n (2). 

Therefore, we have that V = T <2i S(3), where 

ro 1 oooo o ooooooooi-i 
0 0000 20 0 000010000 
0 00000 0 800000100 
0 00 20 00 0 001000000 
0 30 0000 0 000000001 
0 00006 0 000100000 

0 00000 8 000001000 

r(3) _ 0 00 25 00 0 010000000 

1 ~ 1000000000000010 • 

0 000011 0 000010000 
0 00000 0 23 00000100 
0 00 11 00 0 001000000 
30 00000 0 000000010 
0 0000 25 0 00 0 1 0 0000 
0 00000 23 000001000 

L 0 00600 0 0 1000000 oJ 

The required number of multiplications is 10 (4 pre¬ 
multiplications, 2 multiplications in column 4, 2 multiplications 
in column 5, 1 multiplication in column 7 and 1 multiplication in 
column 8). The number of additions is 60 (44 pre-additions, 16 
post-additions). A complexity comparison of A-point FFHT fast 
algorithms (for N = 8 and N = 16) is given in Tables [2] and[3] 


Table 2: Complexity of the 8-point FFHT 


Fast algorithms 

M( 8) 

A(8) 

M(8)+A(8) 

Cooley-Tukey-4 

12 

48 

60 

Split-Radix 

8 

42 

50 

Cooley-Tukey-2 

4 

26 

30 

Rader-Brenner 

2 

24 

26 

Proposed 

2 

22 

24 


Table 3: Complexity of the 16-point FFHT 


Fast algorithms 

M(16) 

A(16) 

Af(16) +A(16) 

Cooley-Tukey-2 

20 

74 

94 

Cooley-Tukey-4 

14 

70 

84 

Split-Radix 

12 

64 

76 

Rader-Brenner 

10 

64 

74 

Proposed 

10 

60 

70 


In the above examples, the Hadamard decomposition 
algorithm presents a lower complexity to compute an 
FFHT compared to existing FFFT/DFT algorithms. Mul¬ 
tiplicative complexity saving regarding classical Cooley- 
Tukey is 50% (N = 16 and N = 8). The total com¬ 
plexity saving regarding the same algorithm is roughly 
25% (N = 16), 20% (N = 8). 

4 Conclusions 

Fast algorithms for the finite field Hartley transform based on 
Walsh-Hadamard decompositions were developed and applied 
to short blocklengths. The theoretical multiplicative complexity 
lower bounds were achieved. The total complexity (additive and 
multiplicative) of the algorithms was compared to that of pop¬ 
ular algorithms and the lower values were obtained for the FT. 
These FTs are attractive and easy to implement using low-cost 
high-speed dedicated circuitry. 
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